## here() starts at /Users/shik544/Projects/hector_project/permafrost_emit
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## Loading hector

1 Introduction

Permafrost—soil that continuously remains below 0°C for at least two consecutive years—underlies an area of 22 (± 3) million km2, roughly 17% of the Earth’s exposed land surface (???), and contains an estimated 1300 (1100-1500) Pg of organic carbon (Hugelius et al., 2014). Recent increases in global air temperature (IPCC, 2013), which are amplified at high latitudes (???), have resulted in widespread permafrost thaw (???), and simulations from variety of climate and land surface models across a wide range of scenarios suggest that this trend will continue into the future (Chadburn et al., 2017; Koven et al., 2013). Permafrost thaw can dramatically alter surface terrain and hydrology (???; ???; ???), with adverse consequences for human infrastructure in permafrost regions (Anisimov and Nelson, 1997; Larsen et al., 2008; Nelson et al., 2002). Moreover, as permafrost thaws, this carbon becomes available to microbes for decomposition, resulting in the production of CO2 and CH4 (???; Brown and Romanovsky, 2008; Romanovsky et al., 2010) that could lead to further warming (Koven et al., 2011; Schuur et al., 2015). Accounting for this permafrost carbon-climate feedback generally increases projections of greenhouse gas concentrations and global temperatures (REFS) and increases estimates of the economic impact of climate change (Chen et al., 2019; Hope and Schaefer, 2015; Yumashev et al., 2019).

The higher carbon emissions associated with warming-driven permafrost thaw may be offset by increases in primary productivity. Some studies based on meteorological tower measurements of carbon flux and optical remote sensing suggest that high-latitude ecosystems (mainly tundra and boreal forest) remain carbon-neutral or are even a small net carbon sink (???; ???). However, other studies, based on different sets of flux measurements and airborne gas sampling, suggest that the high-latitude regions are a net carbon source (???; ???; ???). The uncertainty in the net high-latitude carbon flux may be driven by the heterogeneity of high-latitude landscapes in terms of vegetation cover, soil properties, topography, and many other factors known to affect both the rate of permafrost thaw and the subsequent carbon flux (???; ???; ???; Grant et al., 2019a, 2019b; Turetsky et al., 2002; Wickland et al., 2006). This uncertainty is also present in recent permafrost modeling studies (Burke et al., 2017; Harp et al., 2016; Ito et al., 2016; Qian et al., 2010).

Land surface models are expensive to run, making it challenging to use them for uncertainty quantification and exploration of alternative policy scenarios. Simple climate models are an alternative. (More on simple climate models). (More on permafrost in simple climate models).

Hector (???; Hartin et al., 2015). In this study, we implement a simple representation of permafrost thaw and evaluate its consequences for global climate in Hector.

2 Methods

2.1 Model description: Hector

Hector (Hartin et al., 2015, p.@hartin_2016_ocean,). Simple climate model.

(TODO: More details on terrestrial C cycle in Hector). The default heterotrophic respiration (\(R\)) scheme for a pool \(p\) (detritus or soil) in Hector:

\[ R_{p} = C_p \times f_p \times Q_{10} ^ \frac{T}{10} \]

2.2 Hector permafrost sub-model

The general principle is that permafrost constitutes an additional reserve of soil carbon that, because it is frozen, is inaccessible to microbes. As permafrost thaws and some fraction of this carbon becomes accessible to microbes, it is transferred from the frozen permafrost pool to the standard soil C pool, where it is decomposed by Hector’s standard decomposition routine. Let \(C_{pf}[t]\) be the permafrost C pool and \(C_{s}[t]\) be the soil C pool at time \(t\) (both in units of Pg C), and let \(\Delta C_{pf}[t]\) be the change in the permafrost C pool at time \(t\). The C consequences of permafrost thaw can therefore be represented as:

\[ C_{pf}[t] = C_{pf}[t-1] - \Delta C_{pf}[t] \]

\[ C_{s}[t] = C_{s}[t] + \Delta C_{pf}[t] \]

Let \(C[0]_{pf}\) be the initial size of the permafrost C pool (Pg C) and \(\Phi[t]\) be the fraction of permafrost remaining (in arbitrary area or volume units) at time \(t\). Assuming a uniform permafrost C density, \(\Delta C_{pf}[t]\) can be expressed as:

\[ \Delta C_{pf}[t] = \Delta \Phi[t] C_{pf}[t-1] \]

To a first approximation, \(\Phi[t]\) is a function of air temperature (\(T[t]_{air}\)). Kessler (2015) assume this relationship is linear. However, because the permafrost area fraction is, by definition, bounded by zero (and 11), and because deeper permafrost thaws more slowly than shallow permafrost, we use a negative logistic relationship instead:

\[ \Phi[t] = 1 - \left( 1 + a \exp\left( b * T_{air}[t] \right) \right) ^ c \]

where \(a\), \(b\), and \(c\) are model parameters. The change in frozen fraction at a given time step, \(\Delta \Phi[t]\), is given by:

\[ \Delta \Phi[t] = \max(\Phi[t] - \Phi[t-1], 0) \]

In other words, permafrost thaw is permanent – once it thaws, it does not re-freeze, even if the temperature drops.

For globally-averaged permafrost, \(a = 2.371\), \(b = -0.676\), and \(c = -3.685\) to most closely reproduce the rate of 0.172 K\(^{-1}\) reported by Kessler (2015) over the range of 0.82 to 4 K above the pre-industrial baseline.

Fraction of permafrost thaw as a function of change in global annual mean air temperature since pre-industrial (1750).

Figure 1: Fraction of permafrost thaw as a function of change in global annual mean air temperature since pre-industrial (1750).

2.3 Configuration

Initial permafrost C is set to 1035 Pg C.

3 Results

3.1 Effect of permafrost C

Effect of permafrost C emissions on scenarios.

Figure 2: Effect of permafrost C emissions on scenarios.

4 Discussion

Rate of permafrost C release also depends on soil moisture conditions – drier soils release C much faster (“carbon bomb”) than wetter soils (“carbon fizz”) (Elberling et al., 2013). Moisture will also affect the balance of aerobic (CO2 release) vs. anaerobic (CH4) C release (Turetsky et al., 2002), to the extent that unclear if anaerobic (wet) areas are C sources or sinks (Wickland et al., 2006). Effects of permafrost thaw on soil moisture are a complex hydrological problem – drainage very sensitive to local (micro-)topography (Wickland et al., 2006). So will vegetation cover (Wickland et al., 2006).

Temperature amplification of permafrost carbon feedback (by 2100) 0.02 to 0.36 °C (Burke et al., 2013; Schneider von Deimling et al., 2015, 2012), or 0.1 to 0.8 °C in (MacDougall et al., 2012, 2013), or 10-40% of peak temperature change (Crichton et al., 2016), or 0.2 to 12% (Burke et al., 2017).

Permafrost carbon has greater impact at lower emissions scenarios (Burke et al., 2017; MacDougall et al., 2012, 2013) .

5 Conclusion

6 Acknowledgments

Funded by EPA grant XXX. Cyberinfrastructure support from Pacific Northwest National Laboratory (PNNL).

7 Miscellaneous notes

7.1 Permafrost emissions scenarios

Digitized scenarios from (Schaefer et al., 2011). SiBCASA model predictions of CO2 emissions (permafrost respiration; \(R_{pc}\); note – no methane!) through 2300. These results were digitized using WebPlotDigitizer (https://apps.automeris.io/wpd/), and interpolated to annual resolution (using R stats::spline function).

Digitized scenarios from (Hope and Schaefer, 2015). CO2 and CH4 emissions from SiBCASA model.

(Schuur et al., 2009) – Estimate 0.8 - 1.1 Pg C yr-1.

Back-of-the-envelope estimates from (Zimov et al., 2006): - 500 Gt C in loess that could be completely emitted by 2100 (plus other C sources). - 10-40 g C m-3 day-1 off the bat, slowing down to equilibrium (?) rate of 0.5-5 g C m-3 day-1 for several years. Combine with data on permafrost spatial extent, density, etc. to generate estimates (but can back-calculate from 500 Gt C above?)

7.2 Parameter calibration

We used the BayesianTools R package (Hartig et al., 2019) for all parameter calibration. The outputs of these calibrations are joint posterior distributions of parameters and their covariances, from which we sample for the sensitivity analysis.

For global parameters, we used the following likelihood:

\[ \log(L) = \sum_s Normal(Hector(\beta, Q_{10}, s) | CMIP5(s), \sigma) \]

where \(s\) is one of the four representative carbon pathways (RCPs), \(CMIP5(s)\) are the CMIP5 global mean outputs for the corresponding variables, and \(\sigma\) is the model error (estimated during the fit). We also used the resulting distributions for \(\beta\) and \(Q_{10}\) for the non-permafrost biome in cases 2 and 3. We feel this is appropriate because the CMIP5 models against which these parameters are calibrated do not include permafrost C feedbacks.

For case 2, we calibrated the permafrost-specific \(\beta\) and \(Q_{10}\) against various literature sources, including:

  • Land surface model simulations (Burke et al., 2017; Hope and Schaefer, 2015; Schaefer et al., 2011).
    • Try to calibrate against NPP and soil respiration if possible
  • Literature surveys (Schaefer et al., 2014)
  • Warming experiments (Wickland et al., 2006)
  • TODO: Others?

Some of these are time series, while others are individual estimates at particular points in time. To give them equal weight in the likelihood, we down-weight the time series likelihoods by the number of time steps.

We derived a distribution for the Arctic warming factor from TODO.

TODO: Table and multi-panel figure of input datasets.

For the \(\alpha\) and \(\phi\) parameters in case 3, we looked at the literature on permafrost methane emissions (e.g., Wickland et al., 2006).

7.3 Other notes

Frozen carbon residence time (FCRt) from (Burke et al., 2017):

\[ FCRt = FCRt0 * exp(-\Delta T / \Gamma) (for \Delta T > 0.2 °C) \]

  • \(FCRt_0\) (years) reflects the stability of permafrost C (length of time that permafrost C is stable when \(\Delta T = 0\))
  • \(\Gamma\) – decay term (°C); temperature change at which “the number of years taken for all of the old permafrost C to be emitted reduces to 1/e of its initial value”

Other Hector parameters to consider.

Variable INI name Description Value
\(f_{nv}\) f_nppv Fraction of NPP C transferred to vegetation 0.35
\(f_{nd}\) f_nppd Fraction of NPP C transferred to detritus 0.60
\(f_{nd}\) Fraction of NPP C transferred to soil 0.05
\(f_{lv}\) f_lucv Fraction of LUC change flux from vegetation 0.10
\(f_{ld}\) f_lucd Fraction of LUC change flux from detritus 0.01
\(f_{ls}\) Fraction of LUC change flux from soil 0.89
\(f_{ds}\) Fraction of detritus C that goes to soil 0.60
\(f_{rd}\) Fraction of respiration C to detritus 0.25
\(f_{rs}\) Fraction of respriation C to soil 0.02

According to (Hartin et al., 2015), these were selected to be “generally consistent with previous simple earth system models (e.g., Meinshausen et al., 2011)”.

7.3.0.0.1 pagebreak

8 References

Anisimov, O. A. and Nelson, F. E.: Permafrost zonation and climate change in the Northern Hemisphere: Results from transient general circulation models, Climatic Change, 35(2), 241–258, doi:10.1023/a:1005315409698, 1997.

Brown, J. and Romanovsky, V. E.: Report from the International Permafrost Association: State of permafrost in the first decade of the 21st century, Permafrost and Periglacial Processes, 19(2), 255–260, doi:10.1002/ppp.618, 2008.

Burke, E. J., Jones, C. D. and Koven, C. D.: Estimating the permafrost-carbon climate response in the CMIP5 climate models using a simplified approach, Journal of Climate, 26(14), 4897–4909, doi:10.1175/jcli-d-12-00550.1, 2013.

Burke, E. J., Ekici, A., Huang, Y., Chadburn, S. E., Huntingford, C., Ciais, P., Friedlingstein, P., Peng, S. and Krinner, G.: Quantifying uncertainties of permafrost carbon-climate feedbacks, Biogeosciences Discussions, 1–42, doi:10.5194/bg-2016-544, 2017.

Chadburn, S. E., Burke, E. J., Cox, P. M., Friedlingstein, P., Hugelius, G. and Westermann, S.: An observation-based constraint on permafrost loss as a function of global warming, Nature Climate Change, 7(5), 340–344, doi:10.1038/nclimate3262, 2017.

Chen, Y., Liu, A., Zhang, Z., Hope, C. and Crabbe, M. J. C.: Economic losses of carbon emissions from circum-Arctic permafrost regions under RCP-SSP scenarios, Science of The Total Environment, 658, 1064–1068, doi:10.1016/j.scitotenv.2018.12.299, 2019.

Crichton, K. A., Bouttes, N., Roche, D. M., Chappellaz, J. and Krinner, G.: Permafrost carbon as a missing link to explain CO2 changes during the last deglaciation, Nature Geoscience, 9(9), 683–686, doi:10.1038/ngeo2793, 2016.

Elberling, B., Michelsen, A., Schädel, C., Schuur, E. A. G., Christiansen, H. H., Berg, L., Tamstorf, M. P. and Sigsgaard, C.: Long-term CO2 production following permafrost thaw, Nature Climate Change, 3(10), 890–894, doi:10.1038/nclimate1955, 2013.

Grant, R. F., Mekonnen, Z. A. and Riley, W. J.: Modeling climate change impacts on an arctic polygonal tundra: 1. Rates of permafrost thaw depend on changes in vegetation and drainage, Journal of Geophysical Research: Biogeosciences, doi:10.1029/2018jg004644, 2019a.

Grant, R. F., Mekonnen, Z. A., Riley, W. J., Arora, B. and Torn, M. S.: Modeling climate change impacts on an arctic polygonal tundra: 2. Changes in CO2 and CH4 exchange depend on rates of permafrost thaw as affected by changes in vegetation and drainage, Journal of Geophysical Research: Biogeosciences, doi:10.1029/2018jg004645, 2019b.

Harp, D. R., Atchley, A. L., Painter, S. L., Coon, E. T., Wilson, C. J., Romanovsky, V. E. and Rowland, J. C.: Effect of soil property uncertainties on permafrost thaw projections: A calibration-constrained analysis, The Cryosphere, 10(1), 341–358, doi:10.5194/tc-10-341-2016, 2016.

Hartig, F., Minunno, F. and Paul, S.: BayesianTools: General-purpose mcmc and smc samplers and tools for bayesian statistics. [online] Available from: https://CRAN.R-project.org/package=BayesianTools, 2019.

Hartin, C. A., Patel, P., Schwarber, A., Link, R. P. and Bond-Lamberty, B. P.: A simple object-oriented and open-source model for scientific and policy analyses of the global climate system - Hector v1.0, Geoscientific Model Development, 8(4), 939–955, doi:10.5194/gmd-8-939-2015, 2015.

Hartin, C. A., Bond-Lamberty, B., Patel, P. and Mundra, A.: Ocean acidification over the next three centuries using a simple global climate carbon-cycle model: Projections and sensitivities, Biogeosciences, 13(15), 4329–4342, doi:10.5194/bg-13-4329-2016, 2016.

Hope, C. and Schaefer, K.: Economic impacts of carbon dioxide and methane released from thawing permafrost, Nature Climate Change, 6(1), 56–59, doi:10.1038/nclimate2807, 2015.

Hugelius, G., Strauss, J., Zubrzycki, S., Harden, J. W., Schuur, E. A. G., Ping, C.-L., Schirrmeister, L., Grosse, G., Michaelson, G. J., Koven, C. D. and et al.: Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps, Biogeosciences, 11(23), 6573–6593, doi:10.5194/bg-11-6573-2014, 2014.

IPCC: Climate change 2013 - the physical science basis, UN Intergovernmental Panel on Climate Change: Fifth Assessment Report, doi:10.1017/cbo9781107415324, 2013.

Ito, A., Nishina, K. and Noda, H. M.: Impacts of future climate change on the carbon budget of northern high-latitude terrestrial ecosystems: An analysis using ISI-MIP data, Polar Science, 10(3), 346–355, doi:10.1016/j.polar.2015.11.002, 2016.

Koven, C. D., Ringeval, B., Friedlingstein, P., Ciais, P., Cadule, P., Khvorostyanov, D., Krinner, G. and Tarnocai, C.: Permafrost carbon-climate feedbacks accelerate global warming, Proceedings of the National Academy of Sciences, 108(36), 14769–14774, doi:10.1073/pnas.1103910108, 2011.

Koven, C. D., Riley, W. J. and Stern, A.: Analysis of permafrost thermal dynamics and response to climate change in the CMIP5 earth system models, Journal of Climate, 26(6), 1877–1900, doi:10.1175/jcli-d-12-00228.1, 2013.

Larsen, P., Goldsmith, S., Smith, O., Wilson, M., Strzepek, K., Chinowsky, P. and Saylor, B.: Estimating future costs for Alaska public infrastructure at risk from climate change, Global Environmental Change, 18(3), 442–457, doi:10.1016/j.gloenvcha.2008.03.005, 2008.

MacDougall, A. H., Avis, C. A. and Weaver, A. J.: Significant contribution to climate warming from the permafrost carbon feedback, Nature Geoscience, 5(10), 719–721, doi:10.1038/ngeo1573, 2012.

MacDougall, A. H., Eby, M. and Weaver, A. J.: If anthropogenic CO2 emissions cease, will atmospheric CO2 concentration continue to increase?, Journal of Climate, 26(23), 9563–9576, doi:10.1175/jcli-d-12-00751.1, 2013.

Meinshausen, M., Raper, S. C. B. and Wigley, T. M. L.: Emulating coupled atmosphere-ocean and carbon cycle models with a simpler model, magicc6 - part 1: Model description and calibration, Atmospheric Chemistry and Physics, 11(4), 1417–1456, doi:10.5194/acp-11-1417-2011, 2011.

Nelson, F. E., Anisimov, O. A. and Shiklomanov, N. I.: Climate change and hazard zonation in the circum-arctic permafrost regions, Natural Hazards, 26(3), 203–225, doi:10.1023/a:1015612918401, 2002.

Qian, H., Joseph, R. and Zeng, N.: Enhanced terrestrial carbon uptake in the Northern High Latitudes in the 21st century from the Coupled Carbon Cycle Climate Model Intercomparison Project model projections, Global Change Biology, 16(2), 641–656, doi:10.1111/j.1365-2486.2009.01989.x, 2010.

Romanovsky, V. E., Drozdov, D. S., Oberman, N. G., Malkova, G. V., Kholodov, A. L., Marchenko, S. S., Moskalenko, N. G., Sergeev, D. O., Ukraintseva, N. G., Abramov, A. A. and et al.: Thermal state of permafrost in Russia, Permafrost and Periglacial Processes, 21(2), 136–155, doi:10.1002/ppp.683, 2010.

Schaefer, K., Zhang, T., Bruhwiler, L. and Barrett, A. P.: Amount and timing of permafrost carbon release in response to climate warming, Tellus B, 63(2), 165–180, doi:10.1111/j.1600-0889.2011.00527.x, 2011.

Schaefer, K., Lantuit, H., Romanovsky, V. E., Schuur, E. A. G. and Witt, R.: The impact of the permafrost carbon feedback on global climate, Environmental Research Letters, 9(8), 085003, doi:10.1088/1748-9326/9/8/085003, 2014.

Schneider von Deimling, T., Meinshausen, M., Levermann, A., Huber, V., Frieler, K., Lawrence, D. M. and Brovkin, V.: Estimating the near-surface permafrost-carbon feedback on global warming, Biogeosciences, 9(2), 649–665, doi:10.5194/bg-9-649-2012, 2012.

Schneider von Deimling, T., Grosse, G., Strauss, J., Schirrmeister, L., Morgenstern, A., Schaphoff, S., Meinshausen, M. and Boike, J.: Observation-based modelling of permafrost carbon fluxes with accounting for deep carbon deposits and thermokarst activity, Biogeosciences, 12(11), 3469–3488, doi:10.5194/bg-12-3469-2015, 2015.

Schuur, E. A. G., Vogel, J. G., Crummer, K. G., Lee, H., Sickman, J. O. and Osterkamp, T. E.: The effect of permafrost thaw on old carbon release and net carbon exchange from tundra, Nature, 459(7246), 556–559, doi:10.1038/nature08031, 2009.

Schuur, E. A. G., McGuire, A. D., Schädel, C., Grosse, G., Harden, J. W., Hayes, D. J., Hugelius, G., Koven, C. D., Kuhry, P., Lawrence, D. M. and et al.: Climate change and the permafrost carbon feedback, Nature, 520(7546), 171–179, doi:10.1038/nature14338, 2015.

Turetsky, M. R., Wieder, R. and Vitt, D. H.: Boreal peatland C fluxes under varying permafrost regimes, Soil Biology and Biochemistry, 34(7), 907–912, doi:10.1016/s0038-0717(02)00022-6, 2002.

Wickland, K. P., Striegl, R. G., Neff, J. C. and Sachs, T.: Effects of permafrost melting on CO2 and CH4 exchange of a poorly drained black spruce lowland, Journal of Geophysical Research: Biogeosciences, 111(G2), n/a–n/a, doi:10.1029/2005jg000099, 2006.

Yumashev, D., Hope, C., Schaefer, K., Riemann-Campe, K., Iglesias-Suarez, F., Jafarov, E., Burke, E. J., Young, P. J., Elshorbany, Y. and Whiteman, G.: Climate policy implications of nonlinear decline of Arctic land permafrost and other cryosphere elements, Nature Communications, 10(1), doi:10.1038/s41467-019-09863-x, 2019.

Zimov, S. A., Schuur, E. A. G. and Chapin, F. S.: Climate change: Permafrost and the global carbon budget, Science, 312(5780), 1612–1613, doi:10.1126/science.1128908, 2006.

8.0.0.0.1 pagebreak

9 Appendix

**Figure 1**: Input parameter distributions for global Hector.

Figure 3: Figure 1: Input parameter distributions for global Hector.

**Figure 2**: Input parameter distributions for Hector with biomes.

Figure 4: Figure 2: Input parameter distributions for Hector with biomes.

#> [1] NA
#> [1] NA

9.0.1 Colophon

This report was generated on 2020-02-05 14:34:29 using the following computational environment and dependencies:

#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value                       
#>  version  R version 3.6.2 (2019-12-12)
#>  os       macOS Mojave 10.14.6        
#>  system   x86_64, darwin15.6.0        
#>  ui       X11                         
#>  language (EN)                        
#>  collate  en_US.UTF-8                 
#>  ctype    en_US.UTF-8                 
#>  tz       America/New_York            
#>  date     2020-02-05                  
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  ! package                * version    date       lib
#>    assertthat               0.2.1      2019-03-21 [1]
#>    backports                1.1.5      2019-10-02 [1]
#>    base64url                1.4        2018-05-14 [1]
#>    BH                       1.72.0-2   2019-12-16 [1]
#>    bookdown                 0.16       2019-11-22 [1]
#>    callr                    3.4.0      2019-12-09 [1]
#>    cli                      2.0.0      2019-12-09 [1]
#>    colorspace               1.4-1      2019-03-18 [1]
#>    crayon                   1.3.4      2017-09-16 [1]
#>    desc                     1.2.0      2018-05-01 [1]
#>    devtools                 2.2.1      2019-09-24 [1]
#>    digest                   0.6.23     2019-11-23 [1]
#>    dplyr                  * 0.8.3      2019-07-04 [1]
#>    drake                  * 7.8.0      2019-12-02 [1]
#>    ellipsis                 0.3.0      2019-09-20 [1]
#>    evaluate                 0.14       2019-05-28 [1]
#>    fansi                    0.4.0      2018-10-05 [1]
#>    farver                   2.0.1      2019-11-13 [1]
#>    filelock                 1.0.2      2018-10-05 [1]
#>    fs                       1.3.1      2019-05-06 [1]
#>    ggplot2                * 3.2.1      2019-08-10 [1]
#>    glue                     1.3.1      2019-03-12 [1]
#>    gtable                   0.3.0      2019-03-25 [1]
#>  P hector                 * 2.5.0      2020-01-07 [?]
#>    hector.permafrost.emit * 0.0.0.9000 2019-09-13 [1]
#>    here                   * 0.1        2017-05-28 [1]
#>    highr                    0.8        2019-03-20 [1]
#>    htmltools                0.4.0      2019-10-04 [1]
#>    igraph                   1.2.4.2    2019-11-27 [1]
#>    knitr                  * 1.26       2019-11-12 [1]
#>    labeling                 0.3        2014-08-23 [1]
#>    lazyeval                 0.2.2      2019-03-15 [1]
#>    lifecycle                0.1.0      2019-08-01 [1]
#>    magrittr                 1.5        2014-11-22 [1]
#>    memoise                  1.1.0      2017-04-21 [1]
#>    munsell                  0.5.0      2018-06-12 [1]
#>    pillar                   1.4.3      2019-12-20 [1]
#>    pkgbuild                 1.0.6      2019-10-09 [1]
#>    pkgconfig                2.0.3      2019-09-22 [1]
#>    pkgload                  1.0.2      2018-10-29 [1]
#>    png                      0.1-7      2013-12-03 [1]
#>    prettyunits              1.0.2      2015-07-13 [1]
#>    processx                 3.4.1      2019-07-18 [1]
#>    ps                       1.3.0      2018-12-21 [1]
#>    purrr                    0.3.3      2019-10-18 [1]
#>    R6                       2.4.1      2019-11-12 [1]
#>    Rcpp                     1.0.3      2019-11-08 [1]
#>    remotes                  2.1.0      2019-06-24 [1]
#>    rlang                    0.4.2.9000 2019-12-16 [1]
#>    rmarkdown                2.0        2019-12-12 [1]
#>    rprojroot                1.3-2      2018-01-03 [1]
#>    rstudioapi               0.10       2019-03-19 [1]
#>    scales                   1.1.0      2019-11-18 [1]
#>    sessioninfo              1.1.1      2018-11-05 [1]
#>    storr                    1.2.1      2018-10-18 [1]
#>    stringi                  1.4.3      2019-03-12 [1]
#>    stringr                  1.4.0      2019-02-10 [1]
#>    testthat               * 2.3.1      2019-12-01 [1]
#>    tibble                   2.1.3      2019-06-06 [1]
#>    tidyselect               0.2.5      2018-10-11 [1]
#>    txtq                     0.2.0      2019-10-15 [1]
#>    usethis                  1.5.1.9000 2019-12-16 [1]
#>    withr                    2.1.2      2018-03-15 [1]
#>    xfun                     0.11       2019-11-12 [1]
#>    yaml                     2.2.0      2018-07-25 [1]
#>  source                        
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  local                         
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  local                         
#>  local                         
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  Github (r-lib/rlang@ec7c1ed)  
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  Github (r-lib/usethis@b2e894e)
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#>  CRAN (R 3.6.2)                
#> 
#> [1] /Users/shik544/R
#> [2] /Library/Frameworks/R.framework/Versions/3.6/Resources/library
#> 
#>  P ── Loaded and on-disk path mismatch.

The current Git commit details are:

#> Local:    master /Users/shik544/Projects/hector_project/permafrost_emit
#> Head:     [9ab788d] 2020-02-05: More revisions to introduction

  1. Technically, permafrost area could increase in the case of cooling temperatures, and therefore the area fraction could be greater than 1. However, because even the most aggressive climate action scenarios show temperatures that stabilize above year 2000, we assume that permafrost area will never grow more than the starting value.↩︎

  2. Kessler (2015) report this as temperature change from year 2000. 0.8 K is the warming since pre-industrial as estimated by the default Hector configuration.↩︎